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When locahzed spins on different d orbitals prefer different types of antiferromagnetic ordering, 
the Hund's rule coupling creates frustration. Using spin- wave theory we study the case of two such 
orbitals on a square lattice coupled through Hund's rule, such that the first one couples antiferro- 
magnetically (AF) more strongly to its nearest neighbors, while the second couples more strongly 
to its next nearest neighbors. We find that the zero temperature phase diagram has four regions, 
one characterized by the familiar {n,-K) AF order, a second by the columnar {n,0) order, a third by 
a canted order and a fourth region where a quantum-disordered state emerges. We comment on the 
possible relevance of these findings for the case of Fe-pnictide based antiferromagnets. 

PACS numbers: 74.70.-b,75.10. Jm,75.40.Gb,75.30.Ds 
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The FeAs layer of the parent compound of the Fe- 
pnictide superconductorsi below approximately 134 K 
undergoes a spin-density-wave (SDW) ordering with a 
small magnetic moment 0.35/is per Fe atom^. The d- 
orbitals of the Fe atom are occupied by several electrons 
and in the limit where the Hund's rule coupling is large 
compared to the nearest neighbor (NN) and next near- 
est neighbor (NNN) antiferromagnetic couplings, which 
is believed to be the case for these materials, we may ex- 
pect a much larger moment^'^'^'i per Fe atom 2.6/xb). 
This, and other considerations, have fueled a belief that 
the magnetism in these materials may be of itinerant 
type^iSiiS. In the present paper we explore the possibil- 
ity that the origin of this significantly reduced moment 
is a result of frustration introduced by the fact that the 
various Fe d-orbitals prefer different and competing type 
of magnetic ordering. 

First, let us consider a simplified model in order to 
introduce the reader to the problem discussed here and 
in order to overview our main findings. The more realistic 
modeliiii^, treated within the spin-wave approximation, 
will be presented below. The Hamiltonian 



Ti- — Ji Si^i ■ Sj_i -f J2 Si_2 • Sj^2 

<i3> 

— Jh ■ ^i,2, 



(1) 



describes two distinct spin operators Si_i and Sj^2 corre- 
sponding to two different d orbitals of the same i*'' Fe 
atom. The spins S^.i interact antiferromagnetically with 
their NN S^^i while the spins Si_2 interact with their 
next NN Sj^2 (along the diagonal of the square). The 
Hund's rule coupling Jh tends to align the spins on the 
same atom. The origin of qualitatively different spin- 
interactions for two different d orbitals is discussed in 
Ref. [l^. When Jh — 0, the spins Sj.i order in the (tt, tt) 
order indicated by the red-spins in Fig. [ifa), while the 
spins Si, 2 order in the (tt, 0) (or (0,7r)) order indicated 
by the blue-color spins in Fig.[TI^a). In the absence of Jh 



any choice of direction of order for either type of spins 
is equally acceptable, because our model is rotationally 
symmetric. When Jh > the canted state of Fig. [Tfb) 
is obtained as a compromise state between the two ex- 
tremes of Fig. [T] (a), by tilting the orientation of the blue 
spins by an angle cf) toward the red and, the red spins 
toward the orientation of the blue spins by an angle 9. 
The two spins "bend" towards each other due to Hund's 
rule coupling. Through this canting there is some gain 
from the term proportional to Jh and some loss due to 
both types of spin-spin interactions. When Ji is not too 
different from 2J2 (See Fig. (2), the classical ground state 
is the canted state of Fig. [TJb) for any value of Jh- 

In this paper we also study the role of quantum fiuctua- 
tions around the classical ground states within spin- wave 
theory. We find large amplitude quantum spin fluctua- 
tions when Ji is sufficiently close to 2 J2 and near or in the 
canted phase. Further, we find that for sufficiently large 
Jh/ Ji (and Jh / J2) there is a quantum critical point near 
JhI J\ ^ ^h/(2J2) — 4 (taking 5*1 — S2) from where a 
region of a quantum-disordered state begins. We discuss 
the consequences of our findings for the magnetic state 
of the Fe-pnictides and possible future neutron scattering 
experiments to search for the canted and the disordered 
states. The simpler well-known Ji — J2 modeli^ii^ii^ is 
obtained from our model in the limit of very large Jh- 
However, in order to explain the observed reduced mo- 
ment, the Ji — J2 model requires fine tuning of the Ji/ J2 
ratio to a value very close to the quantum critical point; 
on the contrary, the present model has a much broader 
parameter range yielding large amplitude quantum spin 
fluctuations necessary to explain the observed reduced 
moment in the Fe-pnictides. 

The problem to be discussed here is a somewhat sim- 
plified version of the general Hamiltonian derived in 
Ref. lllilT2l and a generalization of the Hamiltonian given 
by Eq. 



(2) 
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FIG. 1: The classical ground state of the Hamiltonian given 
by Eq. |3l (a) Jh ~* 0. (b) the canted state. The red-color 
spins are described by Ti'^' and they prefer the (7r,7r) state, 
while the blue-colored spins are described by TY'-'^' and prefer 
the columnar order. The Jh term rotates the spins towards 
each other by angles 6 and (j). Red circles denote Fe atoms. 
Orange (yellow) circles denote As atoms above (below) the 
plane. 

<ij> «ij>> 

where the index v = 1,2,..., 5 refers to the five Fe d 
orbitals. When we consider each of the parts H'-"-' sep- 
arately, if 72"^^ > 12, the (7r,0) order (blue-colored 
spins of Fig. [Ha)) is stable, otherwise within spin-wave- 
theory, the (7r,7r) antiferromagnetic order (red-colored 
spins in Fig. [^a)) takes over. For some of the Fe d or- 
bitals in the FeAs based materials J^2^ > /2, while 
for other d orbitals this condition is not satisfie d^^'^^ . 
For simplicity, we will present here the case of just two 
orbitals such that the first one, i.e., 1^ = 1, satisfies the 
condition Jj^-* < Ji^V-^ for (tt, tt) order, while the 1^ — 2 
orbital satisfies the condition J^^"^ > J^^V^ for the (7r,0) 
order. 

First, notice that the Hamiltonian at the classical level 
for certain range of the coupling constants has a ground 
state shown in Fig. [T] (b). The blue spins of the up- 
sublattice are canted by an angle cj) and the red spins 
by an angle 9 as shown in the figure. The total energy 
difference from the energy of the state of Fig. [TJa) is 
SE = cos{2e)-^cos{2(t))~JHSiS2sin{e+(t)), where 

ai = 2Sf{j[^'' - 2J^^^), a2 = 251(27^^ - J^^), and 51,2 
are the maximum length of the two classical spins. In 
the interval < 9 < tt/2, < < tt/2, there are the 
following extrema of the energy. First, the following two 
trivial solutions {9, phi) = (§,0), and {9,cj)) = (0, ^), 



each of which is a stable absolute minimum, respectively, 
when Ci - C2 > 2, and C2 - Ci > 2, where = ^^^7^- 
When neither of these conditions for trivial solutions is 
satisfied the stable absolute minimum is given by 



sin2(20) = C2 ,\r/ ' ^1=^(20) = 7^ ^i'^^^^/.). (4) 

The classical phase diagram is shown in Fig. [2ja). No- 
tice that for any value of the Jh there is the canted 
phase with the angles given as in Eq. [5] provided that 
the other couplings ai) and a2 are not very different from 
each other, i.e., when they satisfy the condition discussed 
above. If, however, these couplings are very different in 
magnitude, the global ground state is the one preferred 
by the stronger coupling, i.e., if a2 » oii the (tt, 0) or- 
der is the ground state, and when ai » a2 the (tTjTt) 
state wins. Both transition lines separating the canted 
order from the (tt, tt) order, or from the (tt, O) order, are 
lines of second order critical points. 




FIG. 2: (a): The classical phase diagram. Here, = ^^^f^ 

and ai = 2S'?(j{^^ - 2J^^'), az = 251(27^^-7^). (b): The 
phase diagram as determined by the spin-wave approxima- 
tion. There is a critical value of {(^i, ("|) and for ("1 > (^i there 
is a range of C^2 where the magnetic order is destroyed and 
a quantum-disordered phase emerges labeled by a question- 
mark in the figure. 

In order to study the role of quantum fluctuations, we 
first carry out a local rotation of the spin quantization 
axes along the direction of the classical order, i.e., by 
angles 9 and cj) for spins on sublattice A as follows: Sfi = 

sm{9)Sli - cos{9)Sf^i, and Sf'^ = cos{9)SIt^ + sm{9)Sf^i, 
while the y component remains unchanged, because we 
have assumed that the rotation is in the x — z plane (the 
plane of the drawing). The expressions for the second 
component are obtained from the above by replacing 9 
tt/2 — (p. For the sublattices B, C and D, we can still use 
the above expressions but with the angles (0, 9) replaced 
by {TT + (f),TT + 9), {TT — (f>,—9) and {—(f>,TT — 9) respectively. 

In order to apply the spin-wave approximation^^, we 
express the operators Sf ^^, Sf,^ Sf,^ using the spin de- 
viation operators, i.e., S^ ,^ — Si, - a\^^ai^y, Sf,^ — 

(ajj, + Gi^ty) and Sf,, = i\f%-{al^ - a^,,,), for both 
cases of spin "color" i' = 1,2. By substituting these oper- 
ators in the Hamiltonian given by Eq. [31 and keeping up 
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TABLE I: The factors needed in Eq[6] are given below. The 
notation: = cos , i 
used. 



cos ky and c^y 



cos fcj; cos ky is 









(1,1) 


2Si{2 cos'^ (6) +Cy sin^(6l)) 


-25'i (c, + Cy cos2(6')) 


(1,2) 


4S'i(cos^(6l)c^a - cos(26l)) 


-45'i sin^(6')c^y 


(2,1) 


252 (2 sin^ (<;/>) +cos2 (<;/>) c J 


— 25'2(ca; + Cy sin^(0)) 


(2,2) 


4S'2(cos(2(?!>) + c^j; sin^(0)) 


— 45'2 COS^ {(f>)Cxy 



to quadratic terms in spin-deviation operators we obtain 



+ ^(Vkaj^^^ak,2 + /i-c + Wkak_iaLk,2 + ^-c), (5) 



where 



Fk - 
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5i52sin(0 + 0)/5. 



-JH/2/sl52(l + sin(0 + 0)), 
JH/2v/^(l-sin((^ + ^?)), 



where the coefficients are given in Table H] and ak,j/ 
are the Fourier components of the operators a^.^ which 
are defined over the entire Brillouin Zone of the non- 
magnetically ordered system, i.e., — tt < kx,ky < tt. 

There are terms proportional to (S'f ,^S'J^ — Sf^^Sj^^) 
arising from both Ti^") and the Jh term. These terms in 
the spin-wave approximation lead to linear terms in the 
operators ai^i, and aj ^ and they can be eliminated by 
choosing the angles to be those minimizing the classical 
energy. 

Now, the quadratic Hamiltonian given by Eq. [5]can be 
diagonalized by means of a canonical transformation 



Ctk.i, 



a 



(6) 



where the coefficients should be chosen to preserve the 
canonical commutation relations for the boson opera- 
tors. This requires the following normalization condition 



,(m) 12 



,(m)|2 



1. Due to the above condi- 



E 

tion, the requirement for the canonical transformation 
to transform the Hamiltonian ([5]) in a diagonal form as 
follows 



n = c 



2 

E 



(7) 



implies that the eigenfrequencies LUk^i, and eigenvectors 
aj^ ^ are given from the set of equations I?(u;k,iy)x'^'^^ = 0, 



where the matrix 



V 



^^k 

W^k 
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Wk 




VFk 




W^k 






Vk 




4^ + 



,(8) 



k,l' ^k,l' '"k,2' 

and v''^\. Here, we have taken advantage of the relations 



and the components of the vector x*^'^^ are u 



and (.-^,)* 



vl^l- We find that 



^k,i/ 



where = (77k,i + ^k,2)/2 
(4''^)2,and Ak 



(9) 

1(1) n(2) , 



^tf^^k^'')^kWk. The staggered magnetizations along the 
direction of the rotated local coordinate system (by the 
angles 6 and (f) are given by 



Su 



(10) 



k At=l 



In the following discussion and calculations presented 
in the figures we restrict ourselves to the special case 

where J^^^ = J^^^ = and, thus, Ci = (52 J//)/(25i j}^^) 
f 2") 

and C,2 = (5*1 J//)/(4S'2 J2 )■ In the entire non-magnetic 
BZ, there are two spin-wave frequencies, an "acous- 
tic" branch, i.e., the o^k,- which vanishes in the long- 
wavelength limit and the "optical" branch Wk,-!- which is 
constant in the long- wavelength limit and of high energy. 
The acoustic frequencies are shown in Fig. [3] along the 
kx and ky directions keeping the value of C2 constant at 
C2 = 4 and varying the parameter C,i. For C2 = 4 there 
are two critical values of ^1 , namely, C,^ — 2 and C,^ — & 
which define the region of the canted phase. The spin- 
wave velocities along the two directions for Ci > ^2 ^ 2 
are different as expected. 

Notice that at the critical point where, we enter 
the canted order from the (tt, tt) order, the modes at the 
wave vectors (tt, 0) and (0, tt) (Figs. [3]) become soft. We 
note that in the pure NN antiferromagnet these modes 
have maximum frequency. At the critical point Ci — Ci~ 
i.e., at border between the canted phase and the (tt, 0) 
phase, these two modes have zero frequency. 

In Fig. [4] we present the staggered magnetizations m\ 
and along the direction of order for spin Si = S2 = 
1/2. The various lines are obtained by keeping (,2 fixed 
and varying C,i. Notice that while the magnitude of the 
staggered magnetization along the rotated direction is a 
continuous function across the transition to the canted 
phase, there are singularities in its derivative at ^1 = 
C,f = C2 ± 2. These singularities indicated by the open 
circles and open squares are caused by the singularities 
in the abrupt change in the angles 6 and 4>. 
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FIG. 3: The lowest ("acoustic") spin- wave frequency band 
along the kx (negative part of the x-axis) and ky (positive 
part of the x-axis) directions for various values of (i which 
correspond to the (vr, tt) order, the canted phase and the (tt, 0) 
order. These results were obtained using ^2=4 and (i — 
1,2,3,4,6,8. 

Notice that for large enough (2 and for comparable 
value of Ci the staggered magnetization along the local 
polarization axes becomes very small. For large values of 
^1 the minimum occurs at the boundary between the 
canted and the (tt, tt) phase. There is a quantum critical 
point which is attained when the Hund's rule coupling is 
large compared to both J2^'' and Jj^''. This limit is be- 
lieved to be the case for the Fe-pnictides. Our model 
reduces to the familiar Ji — J2 model in the limit of 
Jh 00, however, as Fig. [4] indicates reaching this limit 
requires unrealistically large values of Jh as compared 
to all other couplings. Notice, that the transition to the 
canted phase from the side of the (tt, 0) order occurs be- 
fore the staggered magnetization m\ becomes small, even 
for large values of Jh ■ We find that the reason for the en- 
hancement of quantum fluctuations near the Ci = C2 — 2 
boundary is that the spin-wave velocity for large Jh de- 
creases as we increase Jh/J\^'^ and this is not the case 
case at the Ci = C2 + 2 boundary. Therefore, there is 
a quantum disordered phase shown by the green area in 
Fig. [DJb) which illustrates the phase diagram as modified 
by quantum fluctuations. 

In neutron diffraction from the FeAs based antiferro- 



magnets, the canted state should produce a peak with 
intensity proportional to cos^ cf) at k = (ti", 0) (or (0,7r)) 
which has been observed^ and a peak with low intensity 
proportional to sin^ </) at k = (0, tt) (or (tt, 0)). Therefore, 
if the canting angle (j) is small, the latter peak might be 
more difficult to resolve, and, this requires further de- 
tailed experimental investigation. In addition, the mag- 
netic unit cell of the FeAs plane of the canted phase 
is the same as the structural unit cell. This is so be- 
cause there is an orthorhombic lattice distortion below 
155 K and, further, the As atoms are above and be- 
low the plane formed by the Fe atoms in a checkerboard 
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FIG. 4: Comparison of staggered magnetizations m| and m\ 
for spin S\ — S2 ~ 1/2 and for various values of ("2 as a 
function of ("1. Notice that there is a critical value of Jh where 
for Ci C2 — 2 none of these three forms of order survives. 



pattern. Therefore, diffraction using polarized neutrons 
might be a simple way to probe this canted phase. The 
spin-wave dispersion, which is probed by inelastic scatter- 
ing experiments^ has no distinctly different features 
from that of the (tt, 0) phase (see Fig. S]). 

The properties of the quantum-disordered state which 
emerges from the destruction of the long-range order can- 
not be investigated by the present spin-wave theory ap- 
proach. As found in RefsHSITllIiilll, the values of j}''^ 

and J2'^-* are comparable, and, therefore, this phase may 
be accessible by altering these parameters experimentally 
using pressure or electron/hole doping. 
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